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1. September 1 ~ October 30: Working on SVD method 

We have used Maximum Entropy Method in microwave image reconstruction, but 
it failed. We also tried Singular Value Decomposition, we didn’t get good simulation 
results. We continue to work on SVD and get some good results. 

1) The Algorithm: 

Let’s consider a set of linear algebraic equations looks like this: 

a u x, + a i 2 x 2 +a,j x j +... + a, N x N =b, 

a 21 x, + a 22 x 2 + a 23 x 3 + ... + a 2N x N = b 2 

a 3 i x i +a 32 x 2 +a 33 x 3 +... + a JN x N =b 2 (1) 

a Ml X l +a M2 X 2 +a M3 X 3 +"- + a MN X N = b M 

Here the N unknowns x } , j - 1, 2, ... , N are related by M equations. The coefficients a y 
with i= 1, 2, ... , M and j =1, 2, ... , N are known numbers, as are the right-hand side 

quantities 6, , ii = 1, 2, ... , M. 

Its matrix form can be written as: 

A • x = b (2) 

For solution of the above matrix equation for an unknown vector x, we can 
calculate the matrix A' 1 which is the matrix inverse of a square matrix A, i.e., 
A • A' 1 = A' 1 • A = I , where I is the identity matrix ( all zeros except for ones on the 
diagonal). 

For inversion a matrix, we can use Gauss-Jordan elimination or LU 
decomposition. But in many cases dealing with sets of equations or matrices that are 
either singular or else numerically very close to singular, they fail to give satisfactory 
results. A powerful set of techniques to solve the problem is known as singular value 
decomposition, or SVD. 

Singular value decomposition (SVD) is based on the following theorem of linear 
algebra: Any M x N matrix A whose number of rows M is greater that or equal to its 
number of columns N, can be written as the product of an M x N column-orthogonal 
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matrix U, an N x N diagonal matrix W with positive or zero elements ( the singular 
values ), and the transpose of an N x N orthogonal matrix V. The various shapes of these 
matrices will be made clearer by the following tableau: 



The matrices U and V are each orthogonal in the sense that their columns are 
orthonormal, 
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X U *U« 

i=l 

- 6 * 

l<k<N 




l<n<N 

(4) 
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i=l 
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II 

1< k <N 




1< n<N 

(5) 



2) Microwave Imaging: 
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In microwave imaging, the scattered fields can be expressed by the induced 
current distribution using a non-linear integral equation here. 

E = -j<y/i j’jj ) G(x,x , )-j(x , )d 3 x' (8) 

where G is the dyadic Green’s function which can be expressed as: 

G(x, x') = (i - grad. (grad,. <p(x, x’) (9) 

e" jkr 

and O(r) = <D(x, x f ) = - — 

4nr 

and r = |x-x1 = [(x-^ f +(y-r]f +(z-<; f]%, we have 

+ JJ r (- k ’ r ' k ^ 3kr + 3 ) [ J„ ■ (x - x’)+ J, • (y - /)+ J, -(z- Z’) ]e -jkc ds'| 

( 10 ) 

with r = +{y-t\f + (z-^) 2 

and r = (x-^)e I +(y-n)e,+(z-0«« 

In order to use SVD, we can express the equation (10) into matrix form. 
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With (1 1), we can program on SVD to re microwave image. 
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3) Computer simulation and results 



Figure 1 Current and receiving plane Figure 2 Elements on the current plane 

Considering a two-dimensional object situated on the plane XOY (Figure 1) , the 
induced current is distributed on this plane. The receiving surface ( the surface where the 
scattered field is calculated) may be a plane parallel to the object plane with distance r. Or 
it may be the surface of a spherical cap with the radius r. 

In our calculation, the object plane can be divided into I x I square elements as 
shown in Figure 2. Similarly, the receiving surface can also be divided into M x N 
elements according to the requirement of our calculation. 

Figure 3 is the computer flow chart for SVD. 

Using SVD method to reconstruct image under several conditions. The attached 
are the simulation results. 
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Figure 3 Image Reconstruction by SVD 
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By the middle of October, we’ve got some good reconstruction results of 
microwave image. Still in some cases, the results are not satisfied. The problem lies to 
the limited data of measurement surface which doesn’t contain enough information for 
image reconstruction. We tried to get more measured points. But this will need large 
continue computer memory when we solved the inverse problem with SVD. And the 
computer resource is limited. We can not do like this. 

In SVD, we want to solve the AJ by the following equation: 

E g -E 0 =AE = B-AJ (12) 

In the program, we use two-dimension array b(n, n) to express the matrix B NxN . If n is 
'too big, the program doesn’t work due to lack of enough memory. So we try to separate 
the array b(n, n) into n arrays of one-dimension . That is bl(l :n). But this method doesn’t 
work because AJ and A E are all related to the whole array b(n, n). And in Fortran, there is 
no pointer which can do the address operations. This method failed. 

2. November 1 ~ November SO: working on Conjugate Gradient Method 
(CGM) 

We have worked on this method before. But it didn’t work well. We modified the 
program and debugged it. Still it can not go to convergence. We found that we separated 
the real part and imaginary part away during our calculation. And CGM is also written for 
the real number calculation. So we started to modify our program to be able to do the 
complex number calculation in the front part. Then we will convert the complex 
calculation into real number calculation for CGM. This work has not finished. 

3. December ~ April: write and defend the thesis 

In this period, the graduate student arrange the computer simulation data, write the 
thesis. Finally, she will defend the thesis. 


P. S. The computer simulation results. ( pp. 7- 36 ) 
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Simulation Results 



1. Images with Field Measured at Planar Surface 

When the scattered field is measured at a plane which is parallel to the object 
plane, the images recovered as shown in Figure 1.1 - Figure 1.2 with two different kind 
of original induced current distribution. 

The original induced current plane is divided into M x M square elements. The 
area for each element is Sj = 0.1m x 0.1m . This plane is simplified as the J plane. 

The receiving plane where the scattered field is measured is divided into N1 x N2 
-square elements. The area for each element is S E =0.1mx0.1m. This plane is simplified 
as the E plane. 

We can see from these graphs that when the initially guessed J is totally different 
from the original current distribution, the images recovered are not good as those whose 
initially guessed J have the similar distribution with the original ones. An important 
factor to this is that when the receiving surface is a plane, the measured data contain not 
enough information to reconstruct the images. 



- 0.4 - 0.4 

Figure 1.1 (a) The original current J with rectangular distribution 

9x9 square elements in J plane 
6x9 square elements in E plane 
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- 0.4 - 0.4 

(b) Guessed J distribution 



(c) Current image with f = 0. 1 GHz, z = 1 X 
Figure 1.1 Continued 


8 



(b) Guessed J distribution 
Figure 1 .2 Continued 
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-OJ -0.3 


(c) Current image with f = 0. 1 GHz, z = 1 X 
With guessed J in (b) 



-0.25 ^ • 

-0.3 -0.3 


(d) Guess J distribution 
Figure 1.2 Continued 
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1.23 ' ' 

-0.3 -0.3 


(e) Current image with f = 0. 1 GHz, z = 1 A. 
With guessed J in (d) 



-0.25 ^ ■ 

-0.3 -0J 


(f) Guessed J distribution 
Figure 1.2 Continued 
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- 0.3 - 0.3 


(g) Current image with f = 0. 1 GHz, z = 1 X 
With guessed J in (f) 

Figure 1.2 Continued 


2. Images with Field Measured at Spherical Surface 

As we can see from the graphs in section l , the images recovered are not ideal. 
We change the receiving surface from the plane into a spherical cap surface with the 
radius is r as shown in Figure 2. 1 With the spherical cap receiving surface, measured 
points with different coordination system parameters (r, 0, <p) contain more useful 
information to recover the images than those of receiving plane. 

With different parameters such as frequency f, 0 0 , and r, and current J 
distributions, we recover the microwave images as shown in the following sections. 
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1) Simulation at Different Frequency 

In these section, the original current J distribution is | J | = 1 at the center element, 
and | J | = 0 with others. The original current plane is still divided into 9x9 elements 
with area S, = 0.1m x 0.1m of each element. The radius r of the sphere cap equals to IX. 
There are 54 measuring points on the receiving surface with the same r, different cp, and 

different 0 ( for spherical coordinate system) . (0 O ~~) 

As you can see, when the frequency is lower than 0.1 1 GHz, the amplitude of the 
images are very low, but you can get the outlook of the object through the images. When 
the f goes high, the amplitudes of the center point and around all increase. You can hardly 
tell what the original object is like comparing with other J distributions ( Figure 2.3). 
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- 0,4 - 0,4 


Figure 2.2(a) The original current J with one-point distribution 

9x9 square elements in J plane 
54 elements in E surface 



- 0.4 - 0.4 


(b) Guessed J distribution 
Figure 2.2 Continued 
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- 0.4 - 0.4 


(c) Current image with f = 0.0 1 GHz, r = 1 X 



- 0.4 - 0.4 


(d) Current image with f = 0.025GHz, r = 1 X 
Figure 2.2 Continued 
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- 0.4 - 0.4 


(e) Current image with f = 0.05GHz, r = 1 X 



- 0.4 - 0.4 


(f) Current image with f = 0.08GHz, r = 1 X 
Figure 2.2 Continued 
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- 0.4 - 0.4 


(g) Current image with f = 0. 1 GHz, r = 1 X 



- 0.4 - 0.4 


(h) Current image with f = 0.1 1GHz, r = IX 

Figure 2.2 Continued 
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- 0.4 - 0.4 


(i) Current image with f = 0. 1 5 GHz, r = 1 X 



- 0.4 - 0.4 


(j) Current image with f = 0.2GHz, r = IX. 
Figure 2.2 Continued 
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2) Simulation with Different 0 Q 

With different 0 O , the receiving elements vary much. When 0 O = y , the 

receiving surface is a hemisphere. As the 0 O becomes small, this surface becomes a 
spherical cap. In this section, we will simulate with different 0 O . The original current 
J distribution is | J | = 1 at the center 9 elements, and | J | = 0 with others. The original 
current plane is still divided into 9x9 elements with area S, = 0.1m x 0.1m of each 
element. The radius r of the sphere equals to IX. There are 54 measuring points on the 
receiving surface with the same r, different cp, and different 0 ( for spherical coordinate 
system) . 

We can see that when 0 O = — , the image is better than anyone else. 



- 0.4 - 0.4 


Figure 2.3 (a) The original current J with rectangular distribution 

With f = 0.1GHz, r = IX 
9x9 square elements in J plane 
54 elements in E surface 
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- 0.4 - 0.4 


(b) Guessed J distribution 



- 0,4 - 0.4 


(c) Current image with 0 0 = n/6 
Figure 2.3 Continued 
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O.S j 



• 0.4 - 0.4 


(d) Current image with 0 O =x/5 



- 0.4 - 0.4 


(e) Current image with 0 O - 7tj 4 
Figure 2.3 Continued 
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- 0.4 - 0.4 


(f) Current image with 6 0 - /r/3 



- 0.4 - 0.4 


(g) Current image with & 0 = 5x/l2 
Figure 2.3 Continued 
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- 0.4 - 0.4 

(h) Current image with 0 O = tt/ 2 
Figure 2.3 Continue 


3) Simulation with Different r 

Different radius r of the sphere cap receiving surface determines which field is 
measured: far field or near field. In this section, we will study the relation of radius r and 
the image recovered. The original current J distribution is | J | = 1 at the central 9 
elements, and | J | = 0 with others. The original current plane is still divided into 9x9 
elements with area S ; = 0.1m x 0.1m of each element. There are 54 measuring points on 
the receiving surface with different <p and different 0 ( for spherical coordinate system) . 

7U 

The radius r of the sphere cap varies and 9 0 - — . 

Increasing r, the amplitudes of the image decrease. When r < \k , the amplitudes 
decrease more than that when r > 1/t . 
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- 0.4 - 0.4 

(c) Current image with r = 0.42, 



- 0.4 - 0.4 


(d) Current image with r = 0.62. 
Figure 2.4 Continued 
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- 0.4 - 0.4 


(e) Current image with r = 0.8/1 



- 0.4 - 0.4 


(e) Current image with r = 1.0A 
Figure 2.4 Continued 
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- 0.4 - 0.4 


(f) Current image with r = 1 2X 



- 0.4 - 0.4 


(g) Current image with r = l AX 
Figure 2.4 Continued 
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(h) Current image with r = 1.6 A 



• 0.4 - 0.4 


(i) Current image with r = 1.8A 
Figure 2.4 Continued 
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4) Simulation with Different Original Induced Current Distributions 

In this section, we recover the microwave image with different original induced 
current J distribution. The original current plane is still divided into 9x9 elements with 


area S, = 0.1m x 0.1m of each element. There are 54 measuring points on the receiving 
surface with different <p and different 0 ( for spherical coordinate system) . The radius r of 


K 


the sphere equals to 1 X with 0 O = — and f = 0. 1GHz . In this section, the guess current 


distribution is the same in all the simulation as Figure 2.5 (b). 

We can see that when the original plane contain more points whose amplitudes 
are not zero, it is easier to reconstruct the image than those have more zero amplitudes. 



- 0.4 - 0.4 


Figure 2.5 (a) The original current J with diagonal distribution 

9x9 square elements in J plane 
54 elements in E surface 
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- 0.4 - 0.4 


(b) Guessed J distribution 



- 0.4 - 0.4 


(c) Current image 
Figure 2.5 Continued 
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- 0.4 - 0.4 


Figure 2.6 (a) The original current J with triangular distribution 

9x9 square elements in J plane 
54 elements in E surface 



- 0.4 - 0.4 


(b) Current image 
Figure 2.6 Continued 
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- 0.4 - 0.4 


Figure 2.7 (a) The original current J with central line distribution 

9x9 square elements in J plane 
54 elements in E surface 



- 0.4 - 0.4 


(b) Current image 
Figure 2.7 Continued 
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- 0.4 - 0.4 


Figure 2.8 (a) The original current J with band diagonal distribution 

9x9 square elements in J plane 
54 elements in E surface 



- 0.4 - 0.4 


(b) Current image 
Figure 2.8 Continued 



- 0.4 - 0.4 


Figure 2.9 (a) The original current J with triagular distribution 

9x9 square elements in J plane 
54 elements in E surface 



- 0.4 - 0.4 


(b) Current image 
Figure 2.9 Continued 
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- 0.4 - 0.4 


Figure 2. 10 (a) The original current J with uniform distribution 

9x9 square elements in J plane 
54 elements in E surface 



- 0.4 - 0.4 


(b) Current image 
Figure 2.10 Continued 
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- 0.4 - 0.4 


Figure 2.1 1 (a) The original current J with rectangular distribution 

9x9 square elements in J plane 
54 elements in E surface 



- 0.4 - 0.4 


(b) Current image 
Figure 2. 11 Continued 
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